Noname manuscript No. 

(will be inserted by the editor) 



Evolution of Non-Equilibrium Profile in Adsorbate Layer 
under Compressive Strain 



Enzo Granato • S.C. Ying 



(N ■ 

o : 
< 

in 



o 

C/3 



I 

C 

O 
o 



> 
in 
oo 
(N 
m 

o 

(N 



X 



Received: date / Accepted: date 

Abstract We investigate the time evolution of an ini- 
tial step profile separating a bare substrate region from 
the rest of the compressively strained adsorbate layer 
near a commensurate to incommensurate transition. 
The rate of profile evolution as a function of the mis- 
match, coverage and the strength of the substrate po- 
tential are determined by Brownian molecular dynam- 
ics simulations. We find that the results are qualita- 
tively similar to those observed for the Pb/Si(lll) sys- 
tem. The anomalously fast time evolution and sharp- 
ness of the non-equilibrium profile can be understood 
through the domain wall creation at the boundary and 
its subsequent diffusion into the interior of the adsor- 
bate layer. 

Keywords Nanotribology • Friction mechanisms • 
Dynamic modeling 
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1 Introduction 

The system of an adsorbate layer with a lattice mis- 
match to the underlying substrate has been extensively 
studied both from the experimental and theoretical side. 
In the monolayer regime, it provides a realization of 
many two dimensional phases pQ and has also been 
used as a simple model for studying sliding friction 
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phenomena on surfaces 013,4,5 . In particular, the rich 
phenomena of commensurate to incommensurate phase 
transition depends on the competition between the ad- 
sorption energy and the strain energy from the lattice 
mismatch. Recently, Tringides et al. [6 , 7 have done an 
extensive study of the Pb/Si(lll) system and found the 
" Devil's Stair Case" within a narrow range of coverage 
as predicted by the theory. For a larger coverage, they 
have further found that the adsorbate layer exhibited 
anomalous mass transport through the time evolution 
of step profile created by desorbing the adatoms in a 
small region. This anomalous transport is character- 
ized by a sharp non-diffusive profile at all times as well 
as a critical coverage below which the rate of the pro- 
file evolution slows down exponentially with decreasing 
coverage. To date, there is still no detailed understand- 
ing of the microscopic origin of this anomalous mass 
transport. Very recently, Huang et al. [8] have used 
molecular dynamics simulations in a generalized one- 
dimensional Frcnkcl-Kontorova model to demonstrate 
that collective diffusion of Pb atoms in the wetting layer 
on Si(lll) via domain wall motion for coverage beyond 
monolayer can lead to fast mass transport. However, 
the microscopic origin of a critical coverage for the pro- 
file evolution is still not understood. In this paper, we 
present a study of a generic two dimensional adsorbate 
system [3 S 4,5 with a well understood transition from 
a commensurate c(2 x 2) phase to an incommensurate 
phase. A non-equilibrium profile is created in the adsor- 
bate layer and the time evolution of the profile is then 
studied numerically via molecular dynamics simulation. 
In particular, we focus on the rate of profile evolution 
as a function of the lattice mismatch, the coverage and 
the strength of the substrate potential. We find that 
the results are qualitatively similar to those observed 
for the Pb/Si(lll) system. The main results can be un- 
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derstood via microscopic mechanisms that are generic 
to an adsorbate layer under a compressive strain near 
a commensurate-incommensurate transition. In partic- 
ular, the time evolution of the non-equilibrium profile 
is governed by the rate of domain wall creation at the 
boundary and its subsequent diffusion into the interior 
of the adsorbate layer. 

2 Model 

We consider a two-dimensional (2D) adsorbate layer 
with particle interactions modeled by the Lennard- Jones 
(LJ) pair potential, 

v(r)=e [(r o /r) 12 ~2(r o /r)% (1) 

with 2 parameters: minimum energy distance r Q and 
amplitude e G . The substrate is represented by a rigid 
sinusoidal periodic potential, 

V(x, y) =V [2- cos(2irx/a) - cos(2iry/a)], (2) 

with amplitude V and minima on the sites of a square 
lattice with lattice spacing a. This model has been used 
previously in the study of nonlinear sliding friction of 
adsorbed overlayers - The system consists of 

5000 — 10000 adatoms in an area of size L x L. Peri- 
odic boundary conditions are applied. Throughout this 
paper, we will choose the amplitude of the substrate 
potential and the substrate lattice spacing as the scales 
for the energy and length and all quantities will be ex- 
pressed in dimensionless form. In Fig. 1, we show the 
ground state energy of the adsorbate system at half 
coverage, 8 = 1/2, for different values of the parame- 
ter r in the LJ potential. For r a < 1.56, the adsor- 
bate layer is in a commensurate c(2 x 2) phase with 
interparticle spacing of y/2. The ground state energy 
is a minimum at tq — 1.45. Changing the parameter 
r amounts to changing the lattice mismatch or chang- 
ing the compressive strain in the commensurate phase. 
For ro > 1.45, the layer is compressed relative to the 
minimum energy configuration but still remains in the 
commensurate c(2x2) phase as seen in the constancy 
of the amplitude of the structure factor at wavevector 
Q = (-7T, 7r). For r > 1.56, the value of the structure 
factor at wavevector Q = (tt, tt) starts to decrease, in- 
dicating the entrance into an incommensurate phase. 

In the next two sections, we will present results 
when the adsorbate atoms are removed from a circular 
region with radius R << L to mimic the experimental 
creation of a step profile [6l[7] via localized desorption. 
We will then follow the temporal evolution of this non- 
equilibrium profile to explore the dynamics of this sys- 
tem. These studies will be performed by varying either 
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Fig. 1 Structure factor peak at wavevector Q = (n, n) and 
ground state energy as a function of the minimum energy 
separation r in the LJ pair potential. For coverage 8 = 0.5. 

the mismatch through the variation of ro, the cover- 
age or the strength of the substrate potential Vo- The 
calculations were performed using Brownian molecu- 
lar dynamics (MD) simulations [3] with the time vari- 
able discretized in units of St = 0.002 — 0.01r m( j, where 
r md = {ma^/Vof/ 2 . 

3 Numerical results and discussion 

3.1 Non-equilibrium profile evolution as a function of 
lattice mismatch 

For studying the temporal evolution of the non-equilibrium 
profile, the adsorbate layer is first equilibrated at finite 
temperature. Then, adatoms are removed from a cir- 
cular region in the center of the system. The radius of 
the hole is small, corresponding to a removal of only 150 
adatoms for the system of 5000 particles. Thus, the den- 
sity of the region outside does not change significantly 
as the non-equilibrium profile approaches equilibrium. 
Configurations are stored at regular time intervals. 

In Fig. 2, we show snapshots for the evolution of 
the non-equilibrium profile for increasing times. The 
parameters are r Q = 1.53, L = 50, temperature T = 0.4 
and coverage 0.5. For this set of parameters, the ini- 
tial equilibrium phase before the creation of the non- 
equilibrium profile region is the commensurate c(2 x 2) 
phase with interparticle spacing We note that the 
profile remains sharp initially as distinct from the typ- 
ical diffuse profile resulting from pure diffusion dynam- 
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Fig. 2 Snapshots for the evolution of the non-equilibrium 
profile for increasing times. For r = 1.53, L = 50, tempera- 
ture T — 0.4 and coverage 9 = 0.5. The labels correspond to 
times in units of 10 4 MD time steps. 



ics. This is qualitatively similar to the experimental ob- 
servation for Pb/Si(lll), although the small radius of 
the hole region prevents us from determining a conclu- 
sive behavior for the longtime region where the noise is 
significant. To study the time evolution of the profile as 
a function of the lattice mismatch, simulations studies 
were also done for different values of ro at half coverage. 
The time r (in MD steps) that it takes to completely 
refill the circular region and return to the equilibrium 
is obtained and the average refilling velocity is given by 
1/t. The number of particles inside the circular region 
and the radius difference (R a — R(t)), where R a is the 
initial radius of the hole, are both increasing functions 
of time. Figs 3 show the behavior of these quantities 
for r Q — 1.51 and 1.56. Note that the time scale for the 
two different values of ro differ by more than an order 
of magnitude showing that the refilling dynamics speed 
up rapidly beyond r a = 1.51. 

In Fig. 4, we plot the average refilling velocity vs 
ro. There is a critical value of ro ~ 1.50 below which 
the refilling velocity (1/t) is too small to measure in 
our simulation study. This value is smaller than the 
critical value of ro ~ 1.55, corresponding to the bulk 
commensurate-incommensurate transition at the tem- 
perature T — 0.4. Thus the initial state of the adsor- 
bate layer is still in the commensurate c(2 x 2) phase 
and the diffusion dynamics in the bulk is exceedingly 
small. This is the first hint that the boundary effects 
play a significant role in the time evolution of the non- 
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Fig. 3 Time dependence of the number of particles inside the 
hole region (R < R a ) and the radius difference R — R(t), for 
T = 0.4, L = 100 and coverage = 0.5. For r = 1.51 (top) 
and r = 1.56 (bottom) 
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Fig. 4 Inverse refilling time 1/t (in units of 10 3 MD time 
steps) as a function of the LJ pair potential parameter r . 
For L = 50, T = 0.4 and coverage 8 = 0.5. 

equilibrium profile. Beyond the value of ro ~ 1.56, the 
refilling velocity rapidly takes off. 

To understand the role of the boundary effect qual- 
itatively, we can look at the known results of a semi- 
infinite one-dimensional Frenkel-Kontorova chain with 
an open boundary at one end [TlfTO] . Under a compres- 
sional strain, the last adatom is pushed out towards the 
next maximum in the substrate potential due to the 
lack of balancing force from the unoccupied site, with 
the neighboring adatoms making small adjustments to 
reach equilibrium. There is an activation barrier for 
the end atom to hop over to the unoccupied side and 
thus creating a "light" domain wall at the boundary. 
For increasing value of the parameter ro, the mismatch 
and the compressive stress increases and the end atom 
gets pushed further out. Correspondingly, the activa- 
tion barrier for the "light" wall decreases and eventu- 
ally this barrier can vanish leading to the spontaneous 
generation of a "light" wall at the boundary. At finite 
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temperatures, the onset of rapid time evolution of the 
non-equilibrium profile should correspond to the value 
of ro at which the barrier is comparable to ksT. Once 
generated, this "light" wall can now move rapidly fur- 
ther into the chain under the mechanism of biased dif- 
fusion due to the compressive stress. The net effect is 
the motion of the boundary by one lattice spacing into 
the open area. The whole process can then repeat itself 
leading to the advancement of the profile and refilling 
of the initially empty region. Obviously, the process will 
stop at some point if the chain is finite. However, for 
a semi-infinite chain or one maintained at a constant 
chemical potential at the far end, the unfilled region can 
be filled at a steady rate via this mechanism. The cre- 
ation of the "light" wall at the boundary and its rapid 
propagation into the interior is the key ingredient in 
understanding the anomalous evolution of the nonequi- 
librium profile as shown in our numerical study. Note 
that at the onset of this rapid anomalous transport, 
normal bulk diffusion rate of domain walls and defects 
are negligible in the interior of the adsorbate still in 
the commensurate state. However, the "light" wall cre- 
ated at the boundary can diffuse inwards because the 
unbalanced boundary stress essentially provides an in- 
ward driving force leading to a biased diffusion. For the 
same reason, the profile is sharp because the dynamics 
of its evolution is governed by a biased diffusion mech- 
anism rather than the normal diffusive behavior as ob- 
served in typical Boltzman-Matano systems |TT]. The 
activation barrier for the "light" wall has been stud- 
ied quantitatively in the ID Frenkel-Kontorova model 
in the continuum limit [TlllOj but its detailed depen- 
dence on the mismatch for a fully discrete 2D model as 
presented in this study is yet unknown. 

We have also studied the time evolution dynamics as 
a function of the strength of the substrate potential as 
shown in Fig. 5. As the substrate potential gets weaker, 
the rate of time evolution of the profile increases. This 
can be interpreted as due to two factors. The first is that 
the activation barrier of the " light" wall at the bound- 
ary decreases with decreasing strength of the substrate 
potential. The second factor is that the biased diffu- 
sion rate of the " light" wall towards the interior is also 
enhanced through the reduction of the substrate poten- 
tial. 



3.2 Non-equilibrium profile evolution as a function of 
coverage 

While the variation of the LJ potential parameter r$ 
and hence the lattice mismatch nicely illustrates the 
role that the boundary effects play in determining the 
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Fig. 5 Inverse refiling time 1/r (in units of 10 4 MD time 
steps) as a function of the strength of the substrate potential 
V . For L = 50, r a = 1.53, T = 0.4 and coverage 9 = 0.5. 
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Fig. 6 Snapshots of the layer at different times during refill- 
ing of the initial hole. For L = 100, r = 1.53, T = 0.4 and 
coverage 9 = 0.508. The labels correspond to times in units 
of 3 x 10 4 MD time steps. 

time evolution of the non-equilibrium profile, experi- 
mentally variation of the lattice mismatch is hard to 
realize. Instead, the easy parameter to vary in a real 
system such as Pb/Si(lll) is the coverage 017] • In this 
section, we study the effect of varying the coverage be- 
yond the ideal commensurate coverage of 9 = 1/2 on 
the refilling dynamics. Shown in the Figs. 6 and 7 are 
the results of the refilling dynamics as a function of 
the coverage for a fixed value of potential parameter 
r = 1.53. 

As shown in Fig. 7, there is a critical coverage of 
9 ps 0.5 at T=0.4, below which the refilling dynamics 
is too slow to measure in our study. Beyond a value of 
9 w 0.54 the refilling dynamics rapidly take off. Quali- 
tatively, the increase of coverage beyond 9 = 1/2 again 
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Fig. 7 Inverse refiling time 1/r (in units of 10 4 MD time 
steps) as a function of coverage 9 for L — 100, r = 1.53 at 
T = 0.4. 



increases the compressive strain and facilitates the cre- 
ation of a light wall at the boundary, as discussed in the 
last section. However, in this case another factor come 
into consideration in determining the refilling dynam- 
ics. Beyond the half coverage 8 = 1/2, the configuration 
in the bulk is now incommensurate and there are either 
defects or domain walls (heavy type) to accomodate the 
extra coverage. For a system below the Aubry transi- 
tion [1211131111] . these defects are pinned and immobile. 
In this phase, the diffusion of the "light" wall created at 
the boundary into the interior will be blocked by these 
defects and the refilling process will stop. As the cov- 
erage increases, the critical potential strength for the 
Aubry transition decreases [12] and these defects and 
domain walls become mobile. At this point, the refill- 
ing dynamics can proceed smoothly just as discussed in 
the last section where only the value of the mismatch 
is varied. Thus, when the coverage is increased beyond 
the ideal commensurate value of 8 — 1/2, the time evo- 
lution of the profile is determined by the convolution of 
both the activation rate of domain walls at the bound- 
ary and the diffusion rate of defects in the interior of 
the adsorbate layer. We like to point out that our use 
of the term "critical coverage" in the present context 
differ somewhat from what was used in the empirical 
description of the data for Pb/Si(lll)[51[7] where the 
critical coverage was used to describe the value of the 
coverage where the rate crossover from a region with 
an activation barrier increasing with decreasing cover- 
age to another region with weak coverage dependence. 
In particular, the details of the rapid increase of the 
profile evolution rate beyond the critical coverage has 
not been investigated experimentally. 



In this paper we have investigated the time evolution 
of non-equilibrium profile in an adsorbate system near 
a commensurate to incommensurate transition. In the 
commensurate c(2 x 2) phase, the mobility and diffusive 
rate for the bulk system is essentially zero. However if 
the absorbate system is under compressive strain due 
to lattice mismatch with the substrate, then the barrier 
for activation of a "light" domain wall is reduced at the 
boundary separating the adsorbate from the unoccu- 
pied region. For increasing strain energy due to either 
increased lattice mismatch or increase of coverage be- 
yond the ideal coverage for the commensurate c(2 x 2) 
phase, the barrier for this boundary domain wall gen- 
eration will be lowered to the point where it can be 
easily activated thermally. The profile can then evolve 
in time through the repeated activation of the boundary 
domain wall and its biased diffusion into the interior of 
the adsorbate layer. This mechanism is similar to the 
one proposed by Huang et al. [H] . The difference is that 
in our study the coverage dependence of the activation 
barrier of the boundary domain walls comes out nat- 
urally as a consequence of the microscopic model and 
is not invoked as an additional assumption. The profile 
remains sharp during its evolution because it is con- 
tinuously being driven by the unbalanced compressive 
stress at the boundary. When the stress at the boundary 
is increased via increasing the coverage beyond the ideal 
commensurate coverage instead of changing the lattice 
mismatch, an extra ingredient is required besides the 
activation of the boundary domain wall. In this case, 
the interior of the adsorbate layer has intrinsic defects 
or domain walls to accommodate the extra adatoms. 
Depending on the substrate potential and the cover- 
age, these defects can be either pinned or depinned. The 
transition between these two phases correspond to the 
"Aubry transition " well studied in the literature [T2l 
[T3TfT4] . When the domain walls are " pinned" , the prop- 
agation of the boundary generated domain wall into 
the interior will be blocked by the pinned domain walls 
and the time evolution of the non-equilibrium profile 
will also stop. The anomalous fast transport can only 
be achieved when the system is in the unpinned side 
of the incommensurate phase. The transition to this 
side of the " Aubry phase" requires the coverage to ex- 
ceed a critical value, leading to an extra requirement 
for the the rapid propagation of the non-equilibrium 
profile besides the activation of the boundary domain 
wall. This aspect also differs from the generalized one- 
dimensional Frenkel-Kontorova model results of Huang 
et al. [H] where domain walls in the bulk is always mo- 
bile and unpinned beyond the monolayer coverage. Al- 
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though the present study is for a specific model system 
with a relatively simple commensurate- incommensurate 
phase transition, the generic features of the profile evo- 
lution is applicable to a wide class of adsorbate layer 
under compressive strain. The mechanism of anoma- 
lous mass transport via activation and propagation of 
boundary domain walls should play an important role 
in tribological systems. 
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